Известия Томского политехнического университета. 2003. Т. 306, № 5 


Баскина, О. А. Богатиков и др. - М.: Наука, 1985. - 367 
с. 


УДК 550.831.05(571.1) 

СТАТИСТИЧЕСКАЯ ПАРАМЕТРИЗАЦИЯ СИММЕТРИЧНЫХ ГЕОФИЗИЧЕСКИХ ОБЪЕКТОВ 

В.Н. Устинова, И.Г. Устинова 

Томский государственный университет 
Е-шаіІ: и8Ііпоѵа@§§Г.І8и.ги 

В строении и иерархической согласованности дискретных геолого-геофизических объектов выявляется определённая закономерность и 
упорядоченность в их размерах. Она обнаруживается в форме структур, в их пространственном расположении и временном следовании. 
Повторяемость форм проявляется и достаточно легко типизируется в морфологии поверхностей и морфологических сочетаниях геофизи- 
ческих полей. Математическая идентификация типового облика структур эффективно выполняется с использованием автокорреляционного 
анализа и фильтров Винера. 


Геолого-геофизические объекты на любом из изу- 
чаемых уровней организации есть система систем и 
могут рассматриваться как упорядоченное множество 
дискретных элементов [1-3]. Любая геолого-геофи- 
зическая моносистема представима как двух, трёх и 
более компонентная с эмерджентными свойствами. В 
вертикальном разрезе дискретность и вложенность гео- 
лого-геофизических систем обнаруживается в наличии 
разнопорядковых циклитов [4]. На поверхности земных 
оболочек она проявляется в мозаичной, но закономерно 
построенной совокупности геологических объектов 
[5]. Дискретность оболочечных объектов имеет прямо- 
угольную симметрию [6, 7], которая осуществляется 
через систему ортогональных трещин; либо - квази- 
концентрическую [8, 9], связанную со структурами 
центрального типа. Прямоугольные и концентрические 
блоки структурно согласованы и являются составны- 
ми частями дискретно-иерархической [10] блоковой 
системы. 

Вещественно-структурные комплексы осадочных 
нефтегазоносных бассейнов, формирующиеся в дли- 
тельной истории геологического развития, имеют ряд 
устойчивых форм пространственно-морфологического 
проявления, в близких морфологических конфигураци- 
ях обнаруживаются в геофизических полях. Типовые 
морфологические сочетания в нефтяной геофизике 
достаточно надёжно выявляются в палеоповерхностях, 
строящихся по сейсмическим данным, картах геофизи- 
ческих параметров (структурные карты, карты энергий 


отражений и др.). 

Устойчивые морфологические типы рельефа и 
потенциальных полей характеризуются: кольцевым; 
вихревым; спиральным; симметрично-сигмоидным [11] 
(взаимодополняемым по положительным и отрицатель- 
ным формам); двух-, трёх-, четырёх-, семилучевым соч- 
ленением [5] и т.д. основных элементов морфологии. 

Для выявления и истолкования аномалий централь- 
но-зонального типа могут использоваться различные 
методы. Среди них, в силу наличия типических соче- 
таний аномальных проявлений, наиболее эффективны 
вероятностно-статистические методы, методы класси- 
фикации и др. 

В рамках эргодичной и стационарной модели анали- 
зируемого поля (будь то геофизическое поле или рельеф 
поверхности) важные сведения о свойствах аномалий 
можно получить по автокорреляционной функции 
(АКФ), энергетическому спектру, математическому 
ожиданию. Наиболее информативными для оценки 
свойств составляющих потенциальных полей и полей 
сейсмических параметров являются функция автокор- 
реляции и энергетический спектр. Среди параметров 
автокорреляционной функции, характеризующих форму 
и поперечные размеры изучаемых объектов, выделяются 
дисперсия, радиус нулевой корреляции [12] и др. Радиус 
нулевой корреляции даёт представление о скорости 
спада АКФ. 

Для разнотипных геологических и геофизических 
объектов двумерные оценки функции автокорреляции 
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Рис. 1. Аксонометрические изображения для структур второго 
порядка по горизонту ІІ а юго-восточной части Западно- 
Сибирской плиты 


Рис. 2. Функции автокорреляции, построенные вдоль окружнос- 
тей различных радиусов для структурных поверхностей 
верхнеюрских отложений территории исследований 


могут служить критерием обнаружения их на структур- 
ных планах и в геофизических полях [12]. Оценка фун- 
кций автокорреляции для морфологических объектов 
различного типа, по причине трехмерности объектов и 
ограниченности пространственных реализаций, пред- 
ставляет ну^орые трудности. 

Пусть - стационарна процесс с 

д у (т] 2 ким ожиданием ™ * дисперсией 

■ I и функцией автокорреляции 

р . - 

М і' -.2 э 


где 1: - моменты изм^у^уий (либо точки измерений), а 
изучаемый процесс наблюдается на отрезке от О 
у }~~ у /( т 1 — ^-'людений являются величины 

1 — ’ * для которых необходф ^пос- 

троить оценку ^( т ) функции автокорреляции ® на 
~™езке значений [0,Т] аргумента т. Так как наблюдения 
(для структурных поверхностей, геофизических по- 
лей), зачастую, являются функцией двух независимых 
аргументов х,у, при оценке двумерной функции авто- 
корреляции и создании эффективных расчётных схем 
приходится делать некоторые упрощения. 

Из обобщения теоремы Биркгофа-Хинчина [13] для 
эргодических процессов (среднее по совокупности реа- 
лизаций совпадает с математическим ожиданием и про- 
цесс имеет конечную дисперсию) следует возможность 
оценивания корреляционной функции по единственной 


реализации процесса. Вычисление функции автокорре- 
ляции по контуру геометрической фигуры (накладывае- 
мой на площадную реализацию анализируемого поля), 
по форме близкой к натурным объектам (окружности, 
спирали Архимеда, розе Гвидо Гранди и др.), позволяет 
несколько упростить расчётные схемы и уменьшить их 
влияние на результаты оцениващ т р т) 

Если изучаемый ,пппттррр наблюдается на 

отрезке от 0 до Т, где “ ^ 5 г - радиус окружности, 
вдоль которой ведутся наблюдения, окружность радиуса 
г можно задать параметрическими уравнениями 

1^= р со 8 т 

| Ѵ = Р 8ІП Т ' 

V / 

Наблюдения в этом случае становятся функци- 
ей одного независимого переменного 1 и для оценки 
автокорреляционной функции можно воспользоваться 
известной формулой [12] 
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Для типизации форм структур и аномалий потенци- 
альных полей авторы прибегли к расчёту функций авто- 
корреляции по функционально заданным направлениям. 
АКФ, вычисленные для нескольких типов палеоповерх- 
ностей (рис. 1), с использованием окружности, спирали 


37 



Известия Томского политехнического университета. 2003. Т. 306, № 5 


Архимеда и роз Гвидо Гранди и др. показали эффектив- 
ность такого подхода. Типичные формы рельефа явля- 
ются устойчивыми морфологическими проявлениями, 
присутствуют в каждой структуре. Однако в рельефе (в 
морфологии поверхности или в геофизическом поле), 
зачастую, достаточно контрастно проявляется один или 
два преобладающих типа. Универсальность строения 
структур обнаруживается также в том, что типические 
образы имеют центрально-зональную симметрию. 
Симметрия имеет место, как в центральном сочетании 
(«наложении») типовых форм, так и в их центрально- 
периферическом следовании. Например, в лучевых 
типах выявляется нарастание от центра на периферию 
структуры количества лучевых элементов [14]. 

Некоторые оценки автокорреляционных функций, 
рассчитанные по окружностям разных радиусов, 
представлены на рис. 2. Расчёты позволили выявить 
возможность максимального приближения автокорре- 
ляционной функции к типовой форме. 

Структурные карты куполовидных поднятий (к. и.) 
построены в интервале глубин 2. . .3 км, с сечением 0,2 
км, в кружках - геометрические образы типовых форм 
рельефа в плане палеоповерхности - в случае совпаде- 
ния геометрического образа поверхности и геометри- 
ческой фигуры, по которой выполняется расчёт АКФ 
(рис. 2, г, е). Преобладающий тип палеоповерхности 
наглядно проявляется в АКФ (рис. 2, а, б, в). В оценках 
функции автокорреляции обнаруживается централь- 
но-периферическая изменчивость формы поднятия, с 
нарастанием (от свода к обрамлению) изрезанности 
структуры (рис. 2, а, д). 

Полученные при типизации функции автокорреля- 
ции, в свою очередь, несут информацию о пространс- 
твенных свойствах геофизических объектов, могут 
использоваться для фильтрации трёхмерных полей. Для 
случайного геолого-геофизического поля автокорреля- 
ционная функция согласно обратному преобразованию 
Фурье представима интегралом 

+СО +СО 

Р фЬ$] = | | л(м, 5 

-со -со 

где функция ^ ,М ' - спектр мощности с пространс- 
твенными частотами X, ц. Спектр мощности, в свою 
очередь, может быть рассчитан из исходной функции с 
использованием преобразования Фурье, то есть 

ц,м) = 1 

(Р / -со-» 


Сделав замену переменных, х = г со8 ср, у = г 8Іп ср, 
А,=со С 08 \|/, р = со 8Іп \|/, в соответствии с требованиями 
полярной системы координат, где г, со и у, ср- радиу- 
сы, частоты и углы наклона полярных осей, получим 
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к внутр^щ™! интеграл представляет собой 
где - функция Бесселя нулевого 


порядка, согласно преобразованию Ханкеля [15] 

Р (I 

Увеличивая размерность пространства и применяя 
аналогичные рассуждения можно получи^ спектр мощ- 
ности в любом заданном пространстве [16] 
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где I - функция Бесселя порядка т, ® - размер- 
ность координатного пространства равновлияющих на 
исследуемое поле неоднородностей. 

В выражении (1), известном как преобразование 
Винера-Хинчина, верхний предел интегрирования 
бесконечный, однако, при анализе геофизических полей 
мы имеем дело с конечными реаі“УТ“ями. Для полу- 
чения наилучшего приближения ^ интегралом (1) 
в конечных пределах размер реализации Ь необходимо 
выбирать больше интервала существенных изменений 
игральной функции. При задании величины Ь, 
% * будет определена тем точнее, чем быстрее с 
ростом т убывает подынтегральная функция. Для дости- 
жения настоящего требования, с учетом вида осцилли- 
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" = «, 2- ТГ_ 

при - при ' при 

[Т ^ для достоверности расчетов достаточно, 
чтобы ® случайного поля большей размерности 


затухала быстрее. При несоблюдех 'У 7 _ овия (2) появ- 
ляются отрицательные значения % 1 ^ не имеющие 
физическо™ ™" т эла. Очевидно, если в предположении о 
Т "ТТ“іе < - реальуг” ~ — верности, рассчитывать 

ѵ ™ то в этом случае ^ должна удовлетворять 
условию (2), и по спектру Винера может быть получена 
информация о пространственных свойствах объектов 
исследования. Отсюда следует критерий определения 
истинной размерности случайного поля: действительная 
размерность М поля равна максимальной і те в 

(1), при которой вычисленные значения еще 

положительны [17, 18]. Таким образом, форму объекта 
можно оценить по параметру М и радиусу нулевой 
корреляции г 0 [19]. 

Настоящий критерий «знака» информативен в слу- 
чае, когда мы имеем дело с реализацией поля конкретг 
ной размерности, но на геофизические поля оказывают 
влияние объекты с различными пространственными 
свойст у" л, на практике наи^^ее вероятны комбинации 
полей ^ в пространстве ♦ Так, в АКФ проявляет- 
ся преобладающая форма объекта, его латеральная и 
пространственно-зональная изменчивость, пространс- 
твенное строение, которые могут анализироваться, в 
том числе, через оценки спектра. 

Спектр мощности $ к (со) является информативным 
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параметром, характеризующим как энергетический 
состав поля, так и частотный. В силу наличия для 
типовых форм структур собственной формы функции 
автокорреляции, для каждой из них свойственна функ- 
ция Бесселя, наилучшим образом аппроксимирующая 
АКФ при расчёте энергетического спектра Б к (со). 

Оценка пространственного типа структуры (гео- 
физического поля) также может быть выполнена по 
критерию «максимума энергии». Максимальные ампли- 
^тт^чергетического спектра должны соответствовать 
^ ” наилучшим образом аппроксимирующей АКФ 
конкретного типа аномального поля. Для каждой из 
типовых форм (рельефа, геофизического поля) получен 
свой индивидуальный спектр мощности, характери- 
зующийся функцией Бесселя определённого порядка 
т с максимальной интенсивностью значений спектра 
в конкретном диапазоне пространственных частот и 
собственной величиной І 0 (второй нуль автокорреля- 
ционной функции). 

Различие в спектрах мощностей для типовых струк- 
тур геолого-геофизических полей, характеризующих 
морфологические особенности структур центрального 
типа, позволяет использовать величину энергии Е, 
рассчитанную для соответствующей структуры поля 
при заданных параметрах т и частоте со, для оценки 
типа геолого-геофизического объекта. Расчёт энергии 
можно осуществлять в скользящем окне по профилю, 
либо по системе радиальных профилей. По результа- 
там расчётов оценивается полная энергия реализаций 
конкретных аномалий в их контуре. Для повышения 
точности оценок анализируемые профильные, либо 
площадные реализации перед расчётом преобразуются 
в аномальный процесс, путём как минимум десятикрат- 
ного повторения реализаций. 

Типизация в рамках данного подхода АКФ по 
профильным реализациям поля и по площади, ап- 
проксимация ее в частотной области отрезками рядов 
Фурье-Бесселя позволили построить для разделения 
геофизических полей ряд фильтров. В общем случае 
отклик линейного фильтра Т(х,у) на исходную функцию 
и(х,у) равен двухмерной свертке [16]. 

Учитывая радиальную симметрию исходного поля, с 
помощью Фурье-Бесселя образов функции на входе и(г) 
и передаточной характеристики фильтра Ь(г), функция 
на выходе Ф ^ _ и ^ ^ Т] ( р ) 

поскольку 

Ф(Ц) =%Ц) Лщ>, 

щ= ^ +ш*, р =лД 2 + Щ,,ІД 

где * ■ ѵ где - проекции 

со на соответствующие оси координат, ^ - преобразо- 
вание Винера-Хинчина. 

Как показали исследования, оптимальные в сред- 
неквадратическом смысле винеровские фильтры часто 
далеки от фильтров действительно оптимальных для 
подавления аддитивной помехи, а иногда могут даже 


ухудшить картину восстанавливаемого пу^т [20]. 
Объясняется это, во-первых, тем, что и п (АКФ 
полезного сигнала и помехи) не достаточны для опи- 
сания статистических отличий помехи от полезного 
сигнала, во-вторых, среднеквадратический критерий 
неполно описывает действие помехи и, в-третьих, такая 
фильтраго" осу/'эывается на непростой задаче опре- 
деления пс и п прямо из исследуемого поля [17]. В 
связи с этим более надежным критерием оптимальности 
является максимизация отношения энергий сигнала и 
помехи путем эмпирического анализа их энергоемкос- 
ти во всей пе^зе анализируемых частот при разных 
гипотезах об 

Таким образом, учитывая оценк^ шетров 
идеальных фильтров, при гипотезе об = * весовая 
функция иде; I ильтра 

^0д=)2 ХЩ, р 

щ < щ< ц^; 

при ” ^ - 

вая функция идеального низкочастотного фильтра, при 

т.гУ- 1 ЩЦ> 


Р 

[21], <Щ< Щч 


Г ЙЧ2рЙ{ 


О » 


т тьтра, 

Р 1 "III < Щ< Щ ■ г 

кесп- 


ні, нг; щ^ щ ; щ , щ 
где ^ ? - полосы частот 

энергоемкости полезной составляющей, 
1Л " - граничные индексы интервалов частот. 

Для фильтрации в волновом поле конкретного объекта 
(согласно проведённому автокорреляционному анализу) 
функция Бесселя может быть заменена на АКФ типового 
объекта. Качество выполненной фильтрации, с исполь- 
зованием фильтра Винера, определяется среднеквадра- 
тической погрешностью на выходе фильтра, возника- 
ющей за счёт неточного задания автокорреляционной 
функции; коррелируемо стью поискового объекта и 
других составляющих геолого-геофизических объектов, 
геометрические образы которых в той или иной мере 
всегда присутствуют в изучаемом поле; погрешностью 
оценки энергетического спектра и др. Основной труд- 
ностью при выполнении фильтрации является то, что 
исходные реализации имеют конечные размеры. 

Так, разработанные для создания и фильтрации го- 
лографических образов алгоритмы [16-18], могут найти 
применение при анализе и обработке сейсмических 
данных, построении трёхмерных моделей сейсмогео- 
логических полей. 
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